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1. INTRODUCTION 


The Geophysical Fluid Flow Cell (henceforth GFFC) is an 
experiment that was designed to take advantage of the zero gravity 
environment of the Spacelab to effect studies of rotating convection in 
a spherical shell with radial gravity. The radial gravity Is possible 
because In the absence of the Imposed uniform gravitational 
acceleration that dominates the terrestrial laboratory, other forces 
that mimic buoyancy can be substituted. In the GFFC, a hemispherical 
shell of dielectric fluid Is subject to a radial electric field. The 
resulting polarization forces are in the radial direction and directly 
proportional to the density of the dielectric fluid. Thus, a radially 
oriented buoyancy force Is simulated In a manner that Is exactly 
analogous to that occuring In planetary atmospheres and oceans. 
Unfortunately, however, this dielectric buoyancy or dielectric gravity 
Is rather weak for applied fields that don't cause breakdown In the 
fluid. Thus it Is necessary to conduct experiments using It in the 
zero-gravity environment. The uniform gravity of Earth, would, for 
example, completely dominate the motions in the GFFC, and these motions 
would not be good models for the geophysically interesting problems. 

The GFFC was designed to generate a dielectric gravity as large 
as possible with available materials and safety factors taken into 
account. This neccessl tat es a small Inner radius of the hemispherical 
fluid shell, since the electric field goes like I/radlus , and a small 
gap width, since this enhances the field for a given applied electric 
potential. On the other hand, fluid motions of Interest (instability, 
bifurcations, and turbulence or chaos) occur at low effective viscosity 
or high Reynolds number. Thus the gap of the hemispherical annulus 
cannot be made too small. The GFFC was optimized for the study of 
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a particular class of problems that are of great scientific interest 
and that allowed the simplest GFFC configuration In terms of 
engineering, state of the art materials, and flow visualization. 

The problems to be studied In the first flight of the GFFC on 
Spacelab 111 are concerned with the effects of rotation on thermal 
convection driven by a statically unstable thermal heaflng. This 
heating Is either constant on each hemisphere, or can vary in 
latitude. Of principal interest are the wavenumber selection 
mechanisms, the orientation (bandedness) of the convection cells, the 
time dependence of the convection, and the possible generation of any 
differential rotation. The technical advantages of studying these 
problems are: 

1. It is easy to get well Into the supercritical region of 
parameter space. This Includes parameters that are well beyond the 
reach of computer flow models with adequate resolution. 

2. The direct buoyancy driven convective flow (as compared with 
slantwise or barocllnlc ins tabllllty) has a large signature In both 
thermal gradients and velocity. Thus it is relatively easy to observe 
using Schlieren and dye deformation techniques. 

3. The expected experiment time scales, while of order 1 hour, 
are not so long as to be Impractical given the orbit times and power 
restrictions on the Spacelab. 

It Is natural to ask the question, could the GFFC be used to 
study any other global flow problems of fundamental scientific 
Interest. Substantial resources have been Invested in the construction 
of the GFFC, and given that refllghts might be made at relatively low 
cost, it would be of great benefit to the community If the GFFC could 
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be used productively in other areas. This study addresses this 
question. Our rational for evaluating the potential of alternative 
experiments (e.g. other than the direct convection experiments, that 
might themselves benefit from reflight) is based on the following 
constraints : 

1. That the experiments make a contribution to fundamental 
scientific issues. 

2. That the flows, at least for some parameters, be beyond the 
foreseeable access by numerical computation. There is now reasonable 
confidence in computer simulations of large-scale hydrodynamic flows up 
to certain resolution. There are questions of numerical viscosity and 
thermal diffusion at the small scales, as well as problems in some 
models at singularities (like the pole), but many models have been 
tested against both other models and terrestrial laboratory 
experiments. Thus we want to be sure we can learn something 
fundamental from the GFFC that we might not expect to get from 
numerical modelling. This is a fuzzy issue because all numerical 
results are to a certain extent model dependent, whereas the 
experimental conditions are well defined with known representations of 
viscosity and diffusion. Nonetheless, we don*t want to do experiments 
that simply are to be used to verify numerical models! 

3. That the experiments be done with minimal modification to the 
existing GFFC. That is, we don’t want to have to rebuild the apparatus 
in any substantial way. The kinds of changes that might be tolerable, 
and done at very low cost, include things like software adjustments, 
insertion of obstacles in the working fluid gap, changing the fluid to 
one with different properties, etc. Certainly any mechanical changes 
to the spheres or to the optics system would be only justified by 
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an extremely strong (almost sensational) proposal. 

This study Is organized In the following manner. Section 2 
outlines the capabilities of the GFFC as configured for the thermal 
convection experiments. This review Is not meant to be exhaustive as 
the engineering details of the GFFC have been published elsewhere. 
Section 3 addresses the question of the expected meridional and zonally 
symmetric circulations one might expect in the GFFC when heated In 
various statically stable configurations. Then in section 4 we study 
the stability of the computed zonal flow in order to demonstrate that 
baroclinlc instability problems are not feasible in GFFC. Section 5 
extends the results of section 3 to suggest an oceanographic experiment 
that seems to satisfy the above constraints. 

2. CHARACTERISTICS OF THE GFFC 

Figure 1 shows a schematic of the GFFC Instrument (from the GFFC 
crew training handbook, INT012). The working spheres are rotated on a 
turntable. They are heated and cooled by thermoelectric devices 
mounted on the turntable, with overall cooling by avionics air. Flow 
visualization is by back-focus Schlleren and photochroralc dye 
injection. Images are recorded on 16mm film through a complex optical 
path. Essentially hemispherical data dbverage Is Implemented. The 
Schlleren measures the radially averaged horizontal temperature field 
as a function of latitude and longitude. The dye lines are Injected 
radially, span the fluid gap, and successive frames of film record 
their deformations. Data gathering sensitivity Is about .5 degrees 
centigrade for the Schlleren , and between about .01 and 
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Figure 1 
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and Icm/sec for the dye lines. The high dye line figure depends on the 
rotation rate. Since dye Is Injected on the opposite side of the 
sphere from the camera, the sphere must rotate 180^ before a photo can 
be taken. At high flow rates, the dye can deform so much as to be 
impossible to photograph in a short time. Thus, 1 cm/sec velocities 
can only be recorded via dye if the rotation period is short (e.g. 2 
sec.). 

Figure 2 shows a cross-section of the working fluid region. The 
fluid is contained within rigid boundaries. The spherical surfaces are 
conductors, while the radial equatorial boundary is an Insulator. The 
inner radius is Rj, and the outer radius is Rq . The gap-width D is 
just the difference in these two radii. The applied temperatures Tq 
and T^ can be programmed to vary with latitude. Typically the pole 
is hotter than the equator when this is done. The following list 
summarizes the capabilities and currently used fluid parameters of the 
GFFC. 

Ri * inner radius =■ 2.402cm 
Rq “ outer radius “ 3.300cm 
D * gap width = .898cra 

V » kinematic viscosity * .0065cm /sec 

p =* basic fluid density » .761 gm/cm 

a * coefficient of thermal expansion ■ .00134 C"^ 

K » thermal diffusivity » .00090 cm^/sec 

fi =» rotation rate (about .025 to 3.14 rad/sec) 

Ti and Tq ; radial gradient T^-T^ = 0 to 30®* 

latitudinal (0) gradient > 0 to 20® 
max-min temperatures at any point 50® and 20® 

* 10® minimum driving required for adequate Schlleren sensitivity 
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It is useful to go through the calculation for the dielectric 
radial gravity g. There is a potential across the jgap 

V*Ri(RQ/r-l)sin{Dt/D 

where r is the radius and o) is the excitation frequency. This leads to 
a radial electric field 

E-RiRoVsinwt/Dr^ 

so that VE^ - 4(RiRoV/D)^/r® . 

Now the effective dielectric gravity is /2, where e is the 

dielectric constant of the silicone oil (2x10“^^ mks), and the underbar 
denotes a time average. Thus 

gem=2e(RlRoVnii8/D)*/pr® 

It is seen that gem> henceforth simply g, goes up as the applied 
voltage squared, and decreases as r^. The 5*th power dependence is not 
significantly different from a 2 ' nd power dependence with the present 
gap. Figure 3 gives the dielectric gravity, scaled by the square of 
the applied rms voltage. The range of voltage in the GFFC is 40v to 
lOkv. This gives a maximum outer value of 1 m/sec • The inner value 
is a factor of 5 bigger. In the following we will take the average 
maximum value to be 2m/sec . 

The other important parameters for flow problems in the GFFC are 
the Prandtl number P,*»7.0, the aspect ratio 3*”Ri/D>2.67 and some 
measure of the viscosity. The Taylor number is T^^f^D'^/v^. A 
related dimensionless quantity of importance is the Ekman number 

Here f is the Coriolis parameter 2Q. Figure 4 shows the 
Inverse Ekman numbers obtainable in GFFC as a function of the rotation 
period of the apparatus T{ 2 * This period is 1 to 233 seconds. The 
smallest Ekman number is of order .001, or /E of order .03. 
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Finally, some measure of the flov strength Is required* A 
suitable one Is the thermal Rossby number based on the applied lateral 
temperature gradient at the boundaries VT, and assuming geostrophlc and 
hydrostatic balance. This Is 

R - gaVTD/f^L^, 

where L is measure of the horizontal length scale of the motion (or of 
VT). The accessible values of R are given In figures 3 and 6 for two 
values of VT and a range of L* It Is noted that strong quasl- 
geostrophlc flows (R about *3 Is my definition of this range) occur 
typically for small L«2D or so. Thus we conclude that really strong 
non-linear motions In a baroclinic flow will happen only if the motions 
are confined to a channel at mid- or high— la t i t ude . 

Given these numbers it Is possible to use previous stability 
results of Antar and Fowlis (1980) or Glesler and Fowlls (1979) to show 
that baroclinic instability might just be possible in a mid-latitude 
channel in the GFFC at the extreme conditions (Isec period, 30 
degrees). However, it Is by no means clear that a simple basic state 
based on the thermal wind will be accurate at these conditions In the 
Interior of the fluid. For this reason, and because of It Interest for 
the oceanographic case of section 5, we proceed to look at the zonally 
symmetric circulation expected in the GFFC with latitudinal boundary 
heating . 

3. SYMMETRIC CIRCULATIONS 

Because small L flows are the only ones that give big enough R 
for non-linear baroclinic instability (at least potentially), we start 
by considering an Idealllzed mid-latitude channel model. This Is shown 
In figure 7. A standard scaling Is used for the basic equations, and 
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it Is easily shown (c.f. Holton, 1979) that a local cartesian 
coordinate system is accurate enough for our purposes. Hence y 
corresponds to latitude, z to radius, and the notion is assumed 
Independent of x. The upper and lower surfaces have applied 
temperature distributions, and the vertical endwalls can be either 
insulating or conducting (how one would make a thermally conducting 
radial wall in the GFFC with the lOkv across the gap is not at all 
clear to me). The horizontal boundaries are viscous no~sllp, but the 
vertical boundaries are just Impermeable in our simple model. We 
consider motions only on an f-plane for now. 

The goal is to compute the zonally (x) symmetric flow in the slot 
as a function of the external parameters. This is a problem of 
long-standing interest in geophysical fluid dynamics, yet a solution 
for arbitrary forcing is not available. Our approach has been 
suggested in bits and pieces by Barcilon and Pedlosky (1967), Daniels 
(1976), Hlgnett et al. (1981), and Killworth and Manins (1980). 

For small R and E the interior of the box will be geostrophic and 
hydrostatic. There will however be Ekman boundary layers along the two 
horizontal surfaces and these will generate Ekman suction velocities 

w«+/(E/2)* 3u/3y (1) 

at z*0,l respectively where u is the x-invarlant zonal wind. This 
suction will persist unless the thermal stratification gets to be 
extremely big (Barcilon and Pedlosky 1967). Thus the interior vertical 
velocity scale is /E compared with the a priori thermal wind scale. 

The X momentum equation then suggests that the interior meridional 
velocity V will be of order E or R/E whichever is less. Thus, on the 
f-plane, the meridional circulation is very weak , much smaller than 
that given by simple continuity scaling where v*>w*/E for order 1 aspect 
ratios . 


II 
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The interior momentum equations In the small R»E limit reduce 
simply to the thermal wind balance, 

3u/8z — 3T/3y (2) 

Now the general thermal equation 

vTy+wTjj - kV^T (3) 

requires v and w, the meridional and vertical velocities. But since v 
is much less than w (taking R*>/E say) we only need to find w(y). But 
continuity suggests 

W2=-Vy“0 

again because of the smallness of v« Thus w(y) can be simply found 
from (1) and (2). It is, for constant g, 

w’(y)-(Tj|i-T^|o)/(2’.' + e(T|i-Tjo)) (4) 

Combining (4) with (3) gives the non-dimensional non-linear equation 
V^T - ew’(y).Tz(y,z) (5) 

The key parameter here is c (not to be confused with the dielectric 
constant used earlier), 

€-RPr//E (6) 

The physical situation described by this model is shown in figure 
8. The zonal flow u is in thermal wind balance, but the thermal 
distribution itself is determined by an advec t i ve-dl f f ue ^ ve balance. 

The thermal wind generates vertical velocities via the Ekman suction, 
and these can compress the diffusive effects on T Into a thermocllnlc 
layer. The meridional velocity v is so weak as to have little 
influence on the temperature distribution problem. It is clear that 
only if e is small will the interior thermal wind scaling hold up! 

Equation (5) is solved by an iterative relaxation procedure. 
Actually a modified version of (5) is used that Includes horizontal 
thermal advectlon in the Ekman layers, but this process is effective 
only on the lOZ level. Figures 9a, b, and c show the resulting 
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thermal distributions for a conducting right hand boundary and an 
Insulating left hand wall, with a cosine gradient of T along the bottom 
and a uniform top temperature. Case (a) shows the T field when the 
advection terms In (5) are arbitrarily suppressed. Figures 9b and 9c 
show how the downwards vertical velocity suppresses the thermal field 
Into a boundary layer along the lower wall. At e » 3 this layer Is 
already quite thin, meaning that the zonal velocity field u, that 
responds to Ty will not fill much of the box. Note that our extreeme 
GFFC paramters give R»Pr*l so e * 1//E *10. It would appear that the 
thermal wind scaling will break, down, and that the zonal flow will be 
weaker and more structured In the Interior than one would think at 
first, based on using the boundary T gradient In the thermal wind law. 

A similar situation occurs with variable g. The model Is a bit 
more Involved for this case and although the derivation Is not 
presented In detail It follows the spirit of that given above. Figure 
10 shows that the predicted zonal flow, for the same boundary 
conditions as above, decreases substantially as c Increases. 

If one looks closely at equation (5) It Is obvious that 
asymmetries In the thermal boundary conditions will accentuate the 
advectlve term, whereas if advectlve effect will 

be minimized. Unfortunately, flux conditions at boundaries are very 
difficult to achieve In practice. However relaxing the end conditions 
to Insulating, and providing a similar gradient of temperature along 
the top and bottom boundaries does help. Figure 11 shows the thermal 
fields for these conditions. It Is clear that at c > 10 the 
thermocllne layer Is not yet as strong as In the previous case at e ■ 

3. Figure 12 shows the predicted zonal flows for this heating. At 
e*30 the Interior vertical shear has. decreased by a factor of 3 or so 
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from what would be expected by thermal-wind scaling. Figure 13 shows a 
typical w distribution. 

Solutions of this simplified symmetric circulation model suggest 
that the thermal wind scaling used In stability analyses of GFFC type 
flows may be substantially In error If the paramter c Is larger than 
about 1. However, this Is a parameter one would like to be large In 
order to promote Instability (strong driving, low friction, see (11)). 

One full numerical computation was run for the zonally symmetric 
flow in the GFFC under the maximum stress conditions (largest Taylor 
number and maximum applied thermal gradients). We wanted to get some 
Idea of how big a zonal circulation or thermal wind would be 
generated. The numerical model treats 1/r^ gravity along with all the 
boundary conditions in the actual GFFC. It Is described In Hart (1976) 
where results from It were successfully compared with terrestrial 
laboratory experiments. Figure 14 shows the resulting steady 
solution. The R In the figure Is actually a Rayleigh number, not the 
thermal Rossby number defined earlier, but the forcing corresponds to 
20^ from equator to pole, and a stable radial gradient of 5^C. The 
thermal contours are Inverted because the Rayleigh number Is negative 
(stably stratified), so to Interpret the temperature It la necessary to 
multiply the values .shown In 14b by -1. A stable pool develops In the 
equatorial 'stratosphere* so there Is a meridional Jet parallel to the 
axis of rotation (14c). The rest of the meridional circulation Is 
confined primarily to boundary layers. The zonal flow is strongest at 
high and mid-latitudes and nearer the inner sphere where g Is biggest. 
The strength of the zonal wind Is about half that given by a thermal 
wind scaling based simply on the boundary condition. Whether or not 
this zonal wind might be barocllnlcally unstable Is the topic of the 


next section 



Figure 8. Physical structure of ccnvection driven by non-uniform heating 
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Figure 14a . 








32 


4. STABILITY OF THE SYMMETRIC CIRCULATION 

It Is the Intent of this section to address the question of 
barocllnlc stability of the basic zonal flows of the previous section. 
In particular we would like to know if a flow can be sufficiently 
unstable that one might expect strongly non-linear motions will 
develop. Not only are these strongly non-linear waves of fundamental 
interest, but in addition, aqy GFFC experiment studying unstable waves 
must be able to observe them using the current visualization system. 
This latter point means that they must have sufficient amplitude. In 
particular, if the driving is done with boundary differences of order 
10^, as seems to be the a priori case, then the wave radially-averaged 
temperature fluctuations should be at least .3^, or more conservatively 
1® . That is, they should have amplitudes of lOZ of ^Thor.* Without 
doing a full numerical model, we can only approximate the amount of 
super-criticality needed to bring about these amplitudes. However, a 
good estimate of the required super-criticality can be obtained by 
estimating the predicted unstable barocllnlc wave amplitudes using the 
weakly non-linear theory of Drazin (1970, 1972). He derives an 
amplitude equation for unstable quasi-geos trophic disturbances on a 
thermal wind. For steady equilibrated waves, the pL'es«u'*e disturbance, 
relative to the basic state, has amplitude 

Ae - (ci/q2)-^*(R“Rc)- V/Rc (7) 

where is the critical thermal Rossby number required for linear 
instability. The Landau constants c^ are functions of E and the 
aspect ratio. When these are evaluated for typical GFFC numbers, the 
coefficient of the super crl tlcality varies between .13 and .22 
depending on the wavenumber. Now the temperature, within the hydro- 
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static approximation Is of form T-p 2 > Thus the radial averaged 
temperature Is Just the difference In pressures at the top and bottom 
of the gap, and the equilibrium amplitude gives the size of the thermal 
disturbance. We see that the qupercr Icall t y needs to be about .5 or so 
for the waves to be observable. To be strongly non-llnear the 
supercr 1 t Icall ty should be order 1. This should be hept In mind when 
evaluating the meaning of regime diagrams based on linear stability 
theory . 

In considering the linear stability problem It Is first useful to 
attempt a qualitative answer to the question of the effect of the hot 
pole In the GFFC. If one looks at the Integral conditions on 
quasl-geos trophic Instability (Charney and Stern, 1962, Pedlosky, 1964) 
one finds that a necessary condition for Instability Is that 

/ j. 

f /dydz li|>^qy/ ( uq-c)^ ] + 

o c 

2 . ^ 

f dy[i|»2uoz/(uo-c)2] l-/dy[i(/2y^2/(uQ_c)2]o-0 (8) 

t' t' 

where ifi Is the perturbation stream! unct Ion, Uq the zonal flow (a 
function of y and z), and qy the basic state potential vortlclty 
gradient 

qy • $L^/U - Uoyy ~ 3/3z*F» 3uo/ 3z (9) 

In which F Is the rotational Froude number. Now If g Is much larger 
near the lower boundary, the lower boundary Integral will dominate the 
boundary terms (the single Integrals) of (8). Since for a hot pole 
Uqz Is negative the boundary Integrals will be positive In sum. Then 
the necessary condition requires that the potential vortlclty gradient 
qy must be negative In the flow Interior for Instability to be 
possible at all! Looking at (9) Is can be seen that for qy to be 
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negative, that the B-effect must be over-ridden by either the second or 
third terms. The second is the barotropic Instability term and can, in 
some flows like Jupiter's cloud bands, exceed B* However we are 
interested in the baroclinic instability. The third term rarely 
exceeds B in the atmosphere of the Earth, and does not come close in 
the GFFC. Thus the suggestion is that the hot pole may stabilize the 
GFFC to baroclinic waves. To make the equator hot and the pole cold, 
allowing for a 20^ latitudinal difference in temperature, would require 
a massive re-design and re-construction of the instrument. Nonetheless 
it might be that the GFFC flows are so unstable even with a hot 
polethat experiments with non-linear baroclinic waves might be 
possible. 

To look at this aspect of the problem we consider the 
mid-latitude channel problem for quasi -geos trophic instability of the 
zonal flow Uq(z) with To(y»z)=To. We use the Barcllon ( 1964 ) 
model where damping occurs in Ekman layers at the top and bottom of the 
domain. The boundary value problem for the stability of waves of form 

i|»=simry*exp(la(x-ct)»Y(z) 

is 

(Uo-c)(9/3z-F(z)8/9z*Y-(a^+TT^)'P)+’Fqy»0 (10) 

where qy = B / U-9 / 9 z • F ( z ) 9 Uq / 9 z and F(z )»»1/R* 9To/ . This 
interior problem must be solved subject to the top and bottom boundary 
conditions that w-wgj^jQan, where this latter velocity is given by the 
version of (1) that includes the total wave vertical vortlclty instead 
of just 9u/9y. In the equation for F(z), we incorporate the variable 
gravity into the basic state stratification as a factor (1+z/B)^. The 
upper and lower boundary conditions are 

(Uo“C)«9y/9z - f9uo/3z + /E ( a^+x ^ ) Y/ RF ( z ) ia/2 » 0 (11) 

at z*0 and 1 respectively. 
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These equations are solved numerically by a shooting method. The 
computation was checked by using the same basic state as Barcllon 
(1964) and comparing with his analytical results. We actually use a 
basic state Uq and BT^/dz from the numerical simulation on the full 
hemslphere, taking a y-average In mid-latitudes to get Uq . The 
resulting stability diagram for wavenumber 3, which was the most 
unstable as a function of E, Is shown In figure 13. We have taken 0 « 

0 here because the calculation Is simpler. Note that ci*.03 is close 
enough to Cj^«0 to be considered a fair approximation to the neutral 
curve. The region of linear Instability* where one might expect to 
find growing baroclinlc waves In the GFFC is to the left of the curve 
at small values of friction. From our previous discussion concerning 
amplitudes and super c ri t icall ty we see that one would need to have an 
Ekman number of less than .02 to observe the waves* and probably around 
•01 to get near the strongly non-linear regime. This Is a factor of 
three smaller than Is possible In the GFFC. The required low values of 
the Ekman number means one would need a spherical annulus of three 
times the GFFC depth to observe strongly non-linear baroclinlc waves at 
the highest rotation rates used by GFFC. 

Thus we conclude that no GFFC experiments using the current cell 
are likely to exhibit strong baroclinlc waves. There are many 
approximations In the analyses presented above (quasl-geos trophic * y 
Independent basic states* etc.) that might change the quantitative 
details a little bit. But on the basis of this and other studies 
(Fowlls, private communication)* It would be un-rewardlng to try and 
use the GFFC for science problems for which baroclinlc Instability Is 


the central Issue. 
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5» OCEANOGRAPHIC EXPERIMENTS 

Given that eddy generation In a stably stratified flow in the 
GFFC is weak, are there any other kinds of fundamental experiments that 
might be done? There has been much recent progress in understanding 
the wind driven ocean circulation (Rhines and Young 1982) and in 
certain problems associated with forced buoyancy driven motions in the 
wind gyre (Luyten, Stommel, and Pedlosky 1983)* However most of the 
theoretical understanding of the wind-^dr iven circulation comes from 
models that contain a predetermined and usually fixed vertical 
stratification. How is this stratification set up and how is it 
maintained? The history of physical oceanography has followed a course 
of studying the wind-driven, or mechanically forced, circulation 
separately from the thermohall ne, or buoyancy forced, circulation. 

This separation is loosely based on a scale separation for the two 
circulations, the buoyancy one being much deeper and of much longer 
time-scale. The long term goal is to understand the two separately, 
and then attack the coupled problem. Compared with the wind-driven 
circulation, the thermohallne circulation has received little attention 
until recently. A review by Klllworth (1979) points to the many areas 
of buoyancy driven flows of critical importance to oceanography that 
sorely lack understanding. In the last couple of years, however, 
interest in the thermohaline circulation has blossomed and now the 
subject has become one of the really hot topics in ocean science. 

Can the GFFC contribute anything to this problem? We have 
already seen that modelling a strong eddying ocean is probably 
unrealistic. Indeed it would be hard for the GFFC to compete with 
eddy resolving ocean circulation models, partly because of the 
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low Reynold's numbers and lack of strong barocllnlc Instability In the 
GFFC, but also because the measurements of the flows are made at rather 
low resolution (dealing with radial averages of temperature, for 
example). 

There are, as Killworth points out, few, if any, prototypical 
large'-scale rotating and thermally driven flows that are well 
understood. Even the theory for the f-plane zonally symmetric 
thermally forced circulations discussed in section 3 Is In a pretty 
primitive state. The nature of the end-wall boundary layers and their 
stability has not been addressed except superficially. Of course 
f-plane convection can be studied in the terrestrial laboratory 
(Hignett et. al.). So it is anticipated that at least one aspect of 
the thermohallne circulation might be approached terrestrially. But 
terrestrial experiments can never address any questions having to do 
with the global or even large-basin thermohallne circulations because 
in the oceans the B-effect (sphericity) plays a crucial role that can 
never be simulated on Earth. A.nd in fact the meridionally bounded 
3~gyre circulations are probably so different from the f-plane ones 
that the insights gained by completing the f-plane convection studies 
will probably be of little use to large-scale oceanographers. 

Consider for example the simple problem of a box gyre forced 
thermally from above (or below) by a latitudinal temperature gradient. 
That is, take the box of figure 7, add two meridonal walls at x«0,L, 
and let H be a linear function of y (the 3-plane). What will the 
circulation look like? Certainly some of the processes discussed in 
section 4 will still be presentf There will be Ekman suction and 
associated thermal advection tending to produce a thermocline. However 
there are some fundamental and perplexing differences. 
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One might first attempt to describe' the thermally driven gyre 
flow as follows. Since the forcing is y~lnd ep endent , the interior 
should be as well. This then requires complicated eastern and western 
boundary layers, whose strucijure is not known and even hard to 
contemplate. Nonetheless if one knows where the boundary layers are, 
at least one could perhaps make a stretched grid numerical model to 
study the problem. However the situation is probably worse still. 

The continuity equation 
Wz+Vy+Ux “ 0 

gives an estimate for v. If zonality is assumed in the interior, and 
the gyre has scale L in y and x, and D in z, then we find 

V “ wL/D » WEkmanl"/0 

since Ekman suction is presumed to be the dominant vertical velocity 
generator outside strong convective sinking regions. On the other hand 
the Interior dynamics are dominated by the vertical vorticity equation, 
that for ocean scales reduces to Sverdrup balance 

fw^ “ 0v 

Letting vg be the first estimate of v, and Vg be the second, it is 
clear that 

Vg/vg « fL/3 - f/6f 

where 5f is the variation of the Coriolis parameter across the gyre. 
This is Rhines* recirculation index, and for typical gyres varies from 
2 to 5. Thus continuity and Sverdrup are Incompatible with the 
Sverdrup velocity Vg being bigger. This requires an adjustment to 
the continuity equation, namely that u^ be included to make up the 
difference. Thus one suspects that the buoyancy driven flow in this 
very simple gyre model will not be zonal at all. 
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There is a real need for fundamental insights into the 
qualitative nature of buoyancy circulations on the Biplane* What are 
the shapes of the circulations for various degrees of forcing? Where 
Si-’e Lhe boundary layers and what are their functions? Where are the 
sinking (or rising) regions? What is the extent of small scale sinking 
in these regions and how is it coupled to the more or less laminar 
large scale circulation? These are just some of the problems that 
could benefit from experimental guidance* 

Can a computer do the job? State of the art now is a perhaps 10 
layer model with relatively course (20X20 or so) horizontal 
resolution. However the buoyancy circulation spin up time is so long 
(800 years) that few runs can be carried out. These models are not 
eddy resolving in any sense of the word. In studying smooth 
circulations, numerical computations probably exceed the science 
potential of the GFFC because of their much greater information 
gathering. However, it should perhaps be mentioned again that the 
laboratory experiment is well-defined in terms of its small-scale 
truncation, and this may be a significant advantage. There is 
something more thought provoking, too, in trying to understand 
observations of a real fluid motion. Where the GFFC exceeds foreseable 
computational efforts is in problems that have small scale convective 
sinking regions. It is unlikely that such narrow plumes as expected on 
the basis of rotating direct-convection models can be modelled 
numerically along with a large laminar but stratified ocean. 

It is clear that the present configuration of the GFFC is 
adequate for such studies. One would need to place a set of meridional 
barriers into the gap, but this is a relatively inexpensive task. 
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As opposed to the barocllnlc Instability case, the fact that the pole 
Is hot makes no fundamental difference here. One would impose a 
latitudinal thermal gradient on the inner sphere and think of it as the 
ore.an The outer sphere would be constant in T such that the 

polar region of each gyre would be statically unstable. A key issue, 
however, is the data acquisition instrumentation. While it certainly 
can see . convective cells and thermal side-wall boundary layers, it only 
measures the radially-averaged temperature field. Many of the major 
Issues in thermocllne theory revolve around the generation of the basic 
stratif ici tlon in the oceans, so one would really like for this aspect 
of the problem to have detailed thermal sections. Unfortunately, given 
the presence of the lOkv electric fields in the GFFC, it is hard to see 
how in-situ temperature sampling would be possible. Constructing 
additional remote sensing is also impractical. On the other hand, the 
dye lines would give some Information on radial structure, and 
certainly would give the sense of the flow in the gyre, the position of 
asymmetries, internal boundary layers, etc. 

Thus the main conclusion of this paper is that the most promising 
class of geophysical flow studies with the GFFC involving stably 
stratified rotating fluids is that concerned with large scale 
buoyancy-driven circulations. The fundamental science issues are 
largely qualitative and don't depend on precise measurements of 
detailed quantitltes like turbulent fluxes or laminar vorticity 
fluxes. The main question is whether or not the GFFC flow 
visualization instrumentation is adequate to resolve the key science 
issues. This problem can be addressed by further study of prototype 
problems like the one posed in this section. With a little probing 


42 


It should be possible to focus In on some more specific questions so 
that one can more precisely determine whether the optical 
Instrumentation and data gathering capability of the GFFC la up to the 
Casks. It is felt that the current Interest In this class of problem 
In the physical oceanographic community Is a sign of their Importance, 
and suggests Chat further study of possible GFFC contributions is 
warranted . 
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